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Abstract 

We study the spectrum of a one-dimensional Dirac operator pencil, with a coupling 
constant in front of the potential considered as a spectral parameter. Motivated by recent 
investigations of graphene waveguides, we focus on the values of the coupling constant for 
which the kernel of the Dirac operator contains a square integrable function. In physics 
literature such a function is called a confined zero mode. Several results on the asymptotic 
distribution of coupling constants giving rise to zero modes are obtained. In particular, we 
show that the properties of this distribution depend in a subtle way on the sign variation 
and the presence of gaps in the potential. We also observe that variable sign potentials 
may produce complex eigenvalues of the operator pencil. Some examples and numerical 
calculations illustrating these phenomena are presented. 



1 Introduction and main results 

1.1 Statement of the problem 

Consider the system of differential equations 



[V(x)-(X-k)]^-^ = 0, 

ax 



- + [V(x) - (X + k)]ip 2 = 0, 



dx 

on R, where k, A are parameters and V is a potential. Equivalently, one may define a self-adjoint 
operator by 

/V 4- k — V \ 
TV = ( ^ V-k) = ~ ia2V + kas + V ' 

where V = and 02,03 are Pauli matrices. Then (flj) becomes the eigenvalue equation 

da; 



Tytp = \tp, where ip ■ . 

For a given potential V, let us set A = and introduce the ^/-spectrum associated with V: 

= {7 G C : G spec(T 7 y)}. 
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Equivalently, £y is the spectrum of the linear operator pencil 7 1— > To + jV. Our goal is to 
understand the properties of £y, such as symmetries, existence of real and complex (non-real) 
eigenvalues, eigenvalue estimates and asymptotics. Similar problems, as well as some other 
related questions, have been studied in a variety of situations in mathematical literature, see 
e.g. [BiLi] . |GGHK55Vj . [E3], [Hi], [Sch] , 

Whilst the general asymptotic behaviour and estimates in our case are generally in line with 



earlier results (see Theorems 1.4 and 1.6 we should note that our methods allow the widest class 
of potentials), some unexpectedly subtle phenomena occur depending on the properties of V. 
In particular, £y may have totally different structure for single-sign and variable-sign potentials 



(compare Theorems 1.5 and 1.7), as well as for potentials having gaps (i.e., whose support is not 



connected) and for no-gap potentials (see Examples 2.5 and 2.6). Also, variable-sign potentials 



can produce some (or even all) non-real eigenvalues, which have not been studied previously 



(see Theorem 1.7 and Example 2.4) 



In physical literature, this problem appears in the study of electron waveguides in graphene 
(see [HRPI IStDoPoj and references therein). Note that the electron dynamics in graphene 
is governed by the two-dimensional massless Dirac operator, and the one-dimensional system 
([!]) is obtained as a result of the separation of variables: the parameter k corresponds to 
the frequency of the wavefunction ip in the direction parallel to the waveguide. From the 
physical viewpoint, of particular interest are solutions tpi,ip2 G L 2 , which are called confined 
modes. Among them, especially important in the study of conductivity properties of graphene 
are zero modes: L 2 -solutions corresponding to A = 0. (See Section 1.7 for discussion of 



modes corresponding to A ^ 0.) Zero-energy states in graphene have also been studied for 
potentials of other types — see, for instance, [BaTiBrl IBrFrj . It was shown in [HRP] that for 
the potential Vhrp(x) = — l/cosh(x) the solutions of the system can be found explicitly in 
terms of special functions. Moreover, there exists an infinite sequence of coupling constants 7 
such that is an eigenvalue of the operator T 7 y HRP . An attempt to formulate and prove precise 
mathematical statements confirming and generalising the results of [HRPj was the starting point 
of our research. 

1.2 Basic results 

To state precise results we need to make some basic restrictions on the local regularity and 
global decay of the potential V. We shall assume all potentials are real valued and locally L 2 . 
Let Vo denote the class of such potentials which additionally satisfy 

II^IIl^x-i^+i) -> as |x| -> 00; 

roughly, V G Vo if it decays at infinity. In the literature, Vo is sometimes denoted as c$(L 2 ). 

We can define the constant coefficient operator To as a multiplication operator in Fourier 
space. If V G Vo, we show that V is a relatively compact perturbation of To, allowing us to 



define Ty as an unbounded self-adjoint operator on L 2 (see Section 3.1 for more details). The 
same construction can be used for complex-valued potentials (although, of course, the resulting 
operator will no longer be self-adjoint); this allows us to consider T 7 y for any 7 G C. Further 
use of the relative compactness of V leads to the following: 

Theorem 1.1. If V G Vo then £y is a discrete subset ofC 

Remark 1.2. Standard spectrum. The (usual) spectrum of the self-adjoint operator To can be 
computed easily by considering it as a multiplication operator in Fourier space; we get 

spec(To)=R\(-|fc|,|fc|) =:A k , 
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while this spectrum is purely absolutely continuous. When V G Vo we have that Ty is a 
compact perturbation of To, so these operators must have the same essential spectrum (see 
(ReSil section XIII.4]); thus 



spec ess (TV) = spec ess (T ) = A k 



(2) 



The operator Ty may have eigenvalues outside but these must be isolated and of finite 
multiplicity (using the fact that we're dealing with a 1-dimensional problem it is not hard to 
show that these eigenvalues must in fact be simple; a somewhat restricted form of this result is 



given in Lemma 4.1) 



In common with other Dirac operators, Ty possesses a number of elementary symmetries 
which lead to symmetries for the set Sy . In particular, if V G Vo then — Sy = Sy = Sy, while 
Sy is unchanged if we replace k with —k in the definition of To. With this last symmetry in 
mind we shall henceforth assume k > 0; this will enable us to simplify the statement of some 
results. 

To obtain estimates for the distribution of points in Sy we impose extra global decay 
conditions on the potential V. Let Vi denote the class of real valued locally L? potentials 
which satisfy 

\V(x)\ dx < +oo; 



/ 

Jr 



that is, we require V to be integrable. Equivalently, we can define Vi = VoPlL 1 . The class Vi 
is sometimes denoted as £ 1 {L 2 ). 

Firstly we consider the number of points of Sy lying inside the disc {z G C : \z\ < R} of 
radius R > 0. 

Theorem 1.3. Suppose V G Vi. Then 

#(Sy n {z G C : \z\ <R})<C \\V\\ L iR 
for any R > 0, where C is a universal constant (we can take C = 4e/ir). 

This result can be generalised (using a rather different approach) to deal with potentials 



V G Vq which have weaker decay than is required to be L , see Theorem 3.15 



Lower bounds which complement the upper bounds given by Theorem 1.3 can also be 
obtained. Restricting our attention to real points we have the following: 



Theorem 1.4. Suppose V G Vi . Then 



#(Syn[o 5J R]) > 



R 



TT 



V(x) dx 



+ o{R) 



as R — > oo, while the same estimate holds for #(Sy n [— R, 0]) (by symmetry). In particular, 
Sy n R contains infinitely many points if L V{x) dx ^ 0. 



1.3 Single-signed potentials 

In general, the set Sy may contain complex eigenvalues (see Section [2] for some examples of 
explicit potentials which illustrate various possible behaviours for complex points in Sy). Note 
that, even though the operator Ty is self-adjoint (recall that V is real valued), it does not 
follow in general that the corresponding operator pencil should have a purely real spectrum. 
However, if V does not change sign (as in the example considered in [HRP]) all eigenvalues of 
the operator pencil are real: 
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Theorem 1.5. IfVGVo is single-signed then Ey C M. 



By symmetry we can write Ey = {±7„ : n G N} where < 71 < 72 < . . . denotes the 
sequence of positive points in Ey, arranged in order of increasing size. The bound in Theorem 



1.4 can be turned into an asymptotics: 



Theorem 1.6. Suppose V G Vi is single-signed. Then 

R 



#(Eyn[o,i?]) 



7T 



V(x) dx 



+ o(i2) 



as R — > 00. IfVis non-zero we can equivalently write 

TT 



In 



W\ 



L 1 



n + o(n) 



W\ 



TT 



R + o(R) 



as n 



00. 



1.4 Anti-symmetric potentials 



For potentials of variable sign the behaviour of the 7-spectrum may be different, in some cases 
quite drastically so. For anti-symmetric potentials we have the following: 

Theorem 1.7. If V G V is anti-symmetric then Ey n R = 0. 

Note that, the 7-spectrum may still contain an infinite number of complex eigenvalues; see 



Example 2.4 below 



The absence of real points in the 7-spectrum together with Theorem 1.6 shows that the 



lower bound obtained in Theorem 1.4 is quite sharp. 

Remark 1.8. It is easy to see that translating a potential V changes the operator T 7 y to 



something which is unitarily equivalent. In particular, Theorem 1.7 also applies to potentials V 
satisfying the condition V(a + x) = — V(a — x) for some a € K and all The translation 

invariance of our problem will also be used to simplify the presentation of some arguments in 
Section rjl 

1.5 Potentials without gaps 

Let BVq denote the class of compactly supported real valued functions of (totally) bounded 
variation. Clearly BVq C Vj while BVq contains compactly supported piecewise constant 
potentials with a finite number of pieces, as well as compactly supported functions in C . We 
say that a potential V € BVq has no gaps if 



co(sup P (y))ny- i (o) 



0. 



where |5| and co(S) denote the Lebesgue measure and convex hull of a set S C R respectively. 
Theorem 1.9. Suppose V G BVq has no gaps. Then 

R 



#(Eyn[o,i?]) 



TT 



V{x) dx 



+ 0(1) 



as R — > 00. The same estimate holds for #(Ey n [-R, 0]) (by symmetry). 
When L V(x) dx = this result simply states that Ey n K is finite. 
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1.6 Discussion 



Our results give information about the asymptotics of the counting function #(Xy n [0, R\) as 
R — > oo. For any V G Vi the results of Section [L2~| give asymptotic upper and lower bounds of 



C- [ \V(x)\dx and - 



V(x) dx 



(3) 



respectively. Using Theorem 1.3 we can take C = 2e, in which case the upper bound is actually 
uniform for R > 0. For an asymptotic upper bound the constant can be reduced to at least 
C = e in general (see Remark 3.14). Theorem 1.6 shows the constant can be reduced further 
to C = 1 for single-signed potentials; in this case the asymptotic upper and lower bounds agree 
and an asymptotic formula for the points in £y is obtained. For variable-signed potentials 
the quantities f^\V\ = ||V||£i and lf K V differ, leading to differences in the upper and lower 
bounds in ([5]) even if we could take C = 1. For no-gap potentials V G BVq Theorem 1.9 



shows that it is the lower bound that actually gives the leading order term in the asymptotics 
of #(EyH [0,-R]) asii^oo. 

The above results may lead to a hypothesis that, in fact, the lower bound always gives the 
leading order term in the asymptotics of the counting function of the 7-spectrum. However, 
this is not the case; for general (variable-signed) potentials the precise asymptotic behaviour 
of #(Xy n [0,-R]) as R — > 00 appears to depend on V in a rather subtle way. In particular, 
this behaviour appears to be sensitive to 'gaps' in the potential, namely intervals where V = 



appearing between components of supp(y). We refer to Examples 2.5 and 2.6 below, as well 



as to Example [Z4] where a related phenomenon is observed for complex eigenvalues. In fact, it 
is probable that one could construct potentials for which the actual asymptotics lies anywhere 
between the extremes provided by ^ (with C = 1). 



1.7 Non-zero modes 

If we consider the eigenvalues of T 7 y as functions of 7, we can view £y as the set of points at 
which these curves cross 0. One could equally consider crossings at any other point A belonging 
to (— k, k) (the spectral gap of the operator To). This leads to consideration of the set 



{7 G C : A G spec(T 7 y)}. 



With some straightforward modifications, most of our analysis for £y can be carried over to 
£a,v for any A G (— k, k). We now summarise the changes to the main results. 

Theorem |1.1| holds for £a,v- For V G Vrj, £a,v ' s still symmetric under conjugation and 
unchanged if we replace k with —k; however, we cannot expect Xa,v to be symmetric about 
in general (this symmetry generalises to — £a,v = ^-\,v)- Theorem 1.7 does not generalise. 

Theorems 1.3 1.4 1.5 and 1.6 hold for Sa,v with two adjustments; firstly, the constant C 
in Theorem 1.4 may depend on A, and secondly, the results for points in Sa y nl" no longer 



follow "by symmetry" (but can be obtained by similar arguments). 



The latter comment also applies to Theorem 1.9 which otherwise holds for Sa,v i n the case 



that f R V ^ 0. When f K V = we need to impose further conditions on A (to ensure we avoid 
limiting values of the eigenvalues of T 7 y as 7 — > ±00; cf. [GGHKSSVj , Theorem 8.2(i)). 



1.8 Organisation of the paper 

Section [2] is devoted to examples. The main arguments, together with a number of auxiliary 
constructions and results, are collected in Section [3) Theorem 1A is essentially standard; its 
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proof appears in Section 3J. Some key ideas from the Prufer method are introduced in Section 
3.3| In particular, we re-characterise the set £y n 1 in terms of a quantity Ay, which is 



closely related to the Prufer argument (see Proposition 3.7). The asymptotic behaviour of 
Ay is described (Proposition 3.8) and leads directly to Theorem 1.4 Theorem 1.6 follows 
from a related argument, together with the additional monotonicity of Ay for single-signed V 
(as described in Proposition 3.10). An alternative approach based on the Birman-Schwinger 



principle that could be used in the case of single-signed potentials is discussed in Remark 3.3 

A uniform bound on Ay (given in Proposition 3.11) leads through several intermediate 
results to the proof of Theorem |1.3| The derivative of Ay is considered in Sections |3.4| and 



3.5 The justification of the monotonicity result (Proposition 3.10) appears in Section 3.4 



while in Section 3.5 the proof of Theorem 1.9 is reduced to some technical estimates (given in 
Proposition 3.18). 

For the sake of clarity the proofs of the results in Section [3] which require more technical 
arguments are deferred to Section |4j In Section 4A we consider Lemma |3~6] Section 4^2 deals 



with Propositions 3.8 and 3.11 and Section 4.3 is devoted to the proof of Proposition 3.18 
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2 Examples 

2.1 General description 

The main purpose of this Section is to illustrate the results stated above. We restrict our 
attention mostly to piecewise constant potentials with compact support; these allow the easiest 
analysis and already demonstrate the full range of effects. Consider points ao < a\ < • • • < a m 
which partition the real line into m finite intervals Ij = (oj_i,Oj), j = l,...,m, and two 
semi-infinite intervals 2_ = (— oo,ao) and 2+ = (%, +oo). Consider a potential 

V(x) = W(x; [a , . . . , a m ]; {t*, . . . , v m }) := \ V >> ^ j = 1, . . . , m, ^ 

10, iGLU I + , 

with some given real constants Vj. On each interval, we need to solve the equations 



Vifj 2 = + 



(5) 



with V(x) = Vj = const, and then match the solutions to ensure continuity at the points aj. 
The following result is straightforward. 

Lemma 2.1. For a given constant potential V(x) = v such that k / ±yu, the system ((5j) has 
the general solution 

AM M = r (i) ( MV^v ~ k2x ) \ r(2) / Jos_(Vl 2 v -^x) \ 
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lfv = 0, this solution can be equivalently written as 



e kx + 



-kx 



(6) 



In both cases are arbitrary complex constants. 



Remark 2.2. Lemma 2.1 or more precisely equation Q, immediately implies that for any 
compactly-supported (not necessarily piecewise-constant) potential V{x), any eigenfunction ip € 
L 2 (R) satisfies, under the assumption k > 0, 



^i(minsupp(V)) = -02(imnsupp(V)), 
in order to match the L 2 solutions at ±oo. 



?/>i(maxsupp(V)) = — V'2(maxsupp(y)), (7) 



Let us return to the case of a piecewise constant pote ntial Q. The solution on each interval 

Vj and C® = cf \ £ = 1,2. 
ip2(o, m ). Together with continuity 



Ij, j = 1, . . . ,m, can be written down using Lemma 
By Remark 2.2 we have V'i(oo) = ^(oq), and tpi(d 



2.1 



with v 



conditions at each dj, j = 1, ... ,m — 1 this leads to the homogeneous linear system of 2m 

(i) 

equations with respect to 2m unknowns C- , I = 1, 2, j = 1, . . . , m. Denote the determinant 
of the corresponding matrix of coefficients by Dy{^f). As we are looking for a non-trivial .re- 
solution tjj, we have 7 G Sy if and only if Dy{^f) = 0. Thus, in each particular case our problem 
is reduced to constructing Dy{^) and finding its real or complex roots. 

2.2 Calculations, graphs, and further observations 

We visualise the real roots of Dy(~f) by simply plotting its graph for real arguments. In the 
complex case we use the phase plot method (see [WeSe]) in which the value of axgDy{^) = 
DyOy) 

— ilog : ——r is plotted using colours from a periodic scale. The roots of DyH) are singular- 
ity (7) I 

ities of arg Dy(-f) and appear on the phase plot as points at which all of the colours converge. 
The colour scale which we use in all such plots is shown in Figure [T] 




■ 



-n 



n 



Figure 1: Phase plot colour scale for the value of arg Dy (7) • 



In the following examples it is convenient to set 



7 



r 



1. 



Also, we remark that our determinants Dy{^) are defined modulo a real or complex scaling 
constant, which we choose for convenience of presentation. 



Example 2.3 (Illustration of Theorems [L5] and [L6]). Set V x {x) := W{x; [-1,1]; {1}). Then 

2(7cos(27) + sin(27)) 



7-1 



As the potential is single-signed, the spectrum T>y is real, as illustrated in Figures [2] and [3j 
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2(7Cos(27) + sin(2 , y)) 

Figure 2: The graph of D Vl (j) — — — — — — — - — — against real 7 for the potential V\(x) 



V(x;[-1, 1];{1}). 



7 -l 
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n 



2 4 6 



10 12 14 



Figure 3: The phase plot of £^(7) for complex 7. 



For large (positive) values of 7, the solutions of Dy 1 (j) = are approximately those of 
cos(27) = 0, i.e., 



7T 



n + o(n) 



L 1 



matching Theorem 1.6 



Example 2.4 (Illustration of Theorem 1.7). Consider a class of anti-symmetric potentials 



V 2 , g (x) := W(x;[-1 - g/2, -g/2, g/2, g/2 + 1]; {-1,0,1}) parametrised by the gap length 
g > 0. Then, up to a multiplication by a non-zero constant, 



2cosh(5r)(7 2 + 1 - cos(27) + 7sin(27)) + 2 sinh(5)(7 2 cos(27) +73111(7)) 



For any g, the potential Vi^ g is anti-symmetric; hence the spectrum Sy 2 is purely non-real and 
LV 2fl (7) does not have any real roots. This is illustrated in Figure [4] for g = and g = 1- 

It turns out that the behaviour of complex eigenvalues for the potentials V2, g differs sub- 
stantially for zero and non-zero gaps g. By a rather intricate asymptotic analysis of the cor- 
responding transcendental equations we can show that the large eigenvalues with positive real 
parts are asymptotically located on the curves 



Im.7 = ± 



In Re 7 



if 5 = 



(8) 
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-1 1 



t • * » ♦ » ♦ • « »\ »\ ,** »", 
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Figure 4: The graphs of Dy 2 (7) against real 7 for the potentials V2, g (x) := W(x; [—1 
5/2,-0/2, g/2, g/2 + !];{-!, 0,1}) with ,g = (solid black line) and g = 1 (dashed blue line). 



and on the straight lines 



Im 7 = ± arcsinh 



sinh g 



if 5 > 0. 



(9) 



Figures [5] and [6] illustrate this behaviour of complex eigenvalues. For comparison, we also plot 
the corresponding curves ([§]) and ([£]); one can see that the asymptotics is accurate even for the 
low eigenvalues. 

Re(y) 




Figure 5: The phase plot of Dy 2 (7) for complex 7. The solid black curves 
behaviour of the imaginary parts of the eigenvalues. 



14 16 

illustrate the asymptotic 



Example 2.5 (Illustration of the gap effect in one-gap non-zero-integral potentials). Consider 
the one-gap potentials V^^{x) := W(x; [—g — 1, —g, 0, 2]; { — 1,0, 1}) parametrised by the gap 
length g. For each of these potentials, f R Vs : g = 1 and ||V3 )g ||£i = 3. Explicit calculation gives, 
modulo multiplication by a constant, 



r 



[((7 2 + 1) cos(7) — cos(37) + 7sin(37)) cosh(g) + cos(37) + 7sin(37)) sinh(g)] 



The graphs of Dy 39 (7) for real 7 and g = or g = 1 are shown in Figure [7j 
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Figure 6: The phase plot of Dy 2 1 {pf) for complex 7. The solid black lines ([9]), g = 1, illustrate the 
asymptotic behaviour of the imaginary parts of the eigenvalues. 




Figure 7: The graphs of -Dy 3g (7) against real 7 for the potentials V3 iS (x) := [— g 

1, — 0, 0, 2]; {—1, 0, 1}) with g = (solid black line) and g = 1 (dashed blue line). 



We can expect asymptotics of the form 



#(£y 3i9 n[o,i?]) = c 9 - + o(i), 

7T 



as i? — )• 00. For the no-gap potential V3 Theorem 1.9 gives such an asymptotics with Co = 1 = 
Jm^3,i- O n tne hand, £V 31 (7) has three times as many real roots as Zfy 30 (7) (for sufficiently 
large 7). This leads to a constant G\ = 3 = HV^oIIl 1 m the asymptotics for the gap potential 



V31; c.f., the discussion in Section 1.5 



This example is just a partial case of a more complicated phenomenon which we now discuss 
very briefly; a rigorous statement and proof lie beyond the scope of this paper. Consider a general 
(not necessarily piecewise constant) one gap compact potential V{x) such that supp(V) = 
Ii U I2, where I\ and I2 are compact intervals separated by a gap of length g > 0, and assume 
additionally that V{x) does not change sign on either Ij. If the signs of V\i 1 and V\i 2 coincide, 
then the asymptotic behaviour of eigenvalues is given by Theorem |1.5[ If, however, the signs of 
V\i x and V\i 2 are different, then the asymptotic behaviour is given by a complicated formula 
which depends not only upon the gap length g and the values of |/j.V| but also upon the 
rationality or irrationality of the ratio of these two integrals. The rigorous approach to this 
involves an intricate analysis based on the ergodic theorem and is being developed in a work in 
progress. 
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Example 2.6 (Illustration of a twin gap effect). The gap dependence illustrated in the previous 
Example can be made even more dramatic if we consider some special potentials. Introduce the 
symmetric twin gap potentials 

V 4 , g (x) := W{x; [-g -2,-g- 1, -1, 1,0 + 1,0 + 2]; {-1, 0, 1, 0, -1}) 

parametrised by the gap length g. Note that f R V^ g = and ||T4,g||ii = 4 for any g > 0. Figure 
[8]shows the real curves Dy ig {^) for g = 0.5 and g = 1. One can see that there are only two real 
eigenvalues for the former, and an infinite number of real eigenvalues for the latter. In a work in 
progress we show that the asymptotics of #(Sy 4 n [0, R\) as R — > oo changes abruptly between 
R 

0(1) and — ||V4o||li + o(R), depending on whether g < go := arctanh(l/\/2) « 0.8814 or 

7T 

g > go, respectively. 



6r 




Figure 8: The graphs of ZV 4 (7) against real 7 for the potentials V4 >g (x) := W(x; [—g — 2,—g — 
1, — 1, 1, g + 1, g + 2]; {— 1, 0, 1, 0, — 1}) with g = 0.5 (solid black line) and g = 1 (dashed blue line). 



Example 2.7 (Potential from [HRP]). Using a complicated explicit solution involving special 
functions, Hartmann, Robinson and Portnoi found that, for the potential V"hrp(:e) = — l/cosh(x) 
and any k > 0, the positive part of the spectrum £y HRP coincides with the set k — h + N. We 
treat this potential using the Prufer method and plot, for real 7, the quantity cos(Ay HRP (7)); 
by Proposition 3.7 7 G ^v HRP H M if and only if cos(Ay HRP (7)) = (see Section 3.3 for the 
definition of Ay HRP and further details). The curves in Figure g (d rawn for k = 1 and k = 1.5) 
illustrate the result of jHRPj . 



3 Main arguments 

In this section we give the arguments for the main theorems based on a series of more technical 
results; the proofs for the latter will be deferred to Section |4j 



3.1 General 

The unperturbed operator To is an unbounded self-adjoint operator on L 2 whose domain is H , 
the Sobolev space of (C 2 valued) functions on M; that is 

dom(To) = H 1 = {ip G L 2 : € L 2 }. 



11 



cos(Ai/ HRp (y)) 




Figure 9: The graphs of cos(Ay HRP (7)) for k = 1 (solid black line) and k = 1.5 (dashed blue line). 



In fact it is straightforward to check that ||T O -0|| 2 = || V-0|| 2 + k 2 \\ip\\ 2 , so ||T V>|| 

is equivalent 

to H^H^i. It follows that To defines an isomorphism H 1 — > L 2 . 

Next we consider multiplication by an element of Vn. Firstly note that a norm can be defined 
on Vn using the expression 

||V||v = sup 11^11^2(^-1,^+1)- 

This norm makes Vn a Banach space in which is a dense subset. 

Lemma 3.1. Multiplication by a fixed V € Vn defines a compact map H 1 — > L 2 . 

Proof. Initially suppose V G . Choose a bounded interval / with supp(V) C I. We can 
view multiplication by V as a composition ip h-» ip\r >->■ (Vtj))\i \-t Vip where we firstly restrict 
to /, then multiply by V and finally extend by 0. This gives a map H 1 — > L 2 (I) — > L 2 {I) — > L 2 , 
where the last two steps are continuous and the first step is compact (by the Rellich-Kondrachov 
Theorem; see [Ad], for example). 

Since Cq° is dense in Vn and the set of compact maps is closed, it now suffices to show that 
multiplication defines a continuous bilinear map Vn x H 1 — > L 2 . To this end firstly note that 
the Sobolev Embedding Theorem (ibid.) gives 

IIVHU^rc-M+l) < C\\lp\\ H i( x _ 1:X+1 ) 

for some constant C (which is independent of x). Thus 

W^W^ix-Lx+l) < II^IIl 2 (x-1,x-+1) \\i>\\L™(x-l,x+l) < C||^IU 2 (a;-l,a+l) IIV' \\m{x-l,x+l) ■ 

On the other hand 

\mh = lJ R \mh(x-i,x + i)^ 

with a similar expression holding for || • Combining the above then gives 

ii^Hi, < ^ 2 ^i^ni 2(:r -i 1 , +1 )iiV'iil fl ( a; - M+1 )dx 

<^ 2 su P ||Fiii 2(:c _ lj:c+1) f \\n 2 HHx -i, x+ i)^ =c 2 \\v\\ 2 m 2 h i, 

for any V € V and i\) G H 1 . □ 
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Since H 1 = dom(To), Lemma 3.1 is equivalent to the statement that (multiplication by) 
V G Vo is a relatively compact perturbation of To- It follows that the sum Tq + V = Ty defines 
a self-adjoint operator with dom(Ty) = dom(To) = H 1 (see [ReSi]). 

Although we're interested in real valued potentials it is helpful to consider some basic results 
for more general complex valued potentials as well. Let Vo = Vo ®M C denote the complex 
valued version of Vq (in other words, Vq consists of functions of the fo rm U + iW where 



U,W G V ). Note that -yV G V for any V G V and 7 G C. Now Lemma |3_l| clearly extends 
to Vo, so (multiplication by) some V G Vo is still a relatively compact perturbation of TV Thus 
the sum To + V = Ty defines a closed operator with dom(Ty) = dom(To) = H l . Although 
the operator Ty will not be self-adjoint (unless V G Vo), the essential spectrum is still given 
by ([2]), while any spectrum of Ty in C \ Aj, consists of isolated eigenvalues of finite (algebraic) 
multiplicity. 

Proof of Theorem \l.l[ We have 7 G Xy iff (To + jV)ij) = for some non-trivial tp G H . In 
turn this is equivalent to VTo -1 <£> = flip where ip = Tqi^ G I? and \i = — 1/7; in other words, 
— 1/7 should be an eigenvalue of the operator VTo -1 . However VTo -1 is a compact operator 
by Lemma [3~T| so must have discrete spectrum away from 0. The result follows. □ 



It is helpful to have a more symmetrised version of the idea just used in the proof of Theorem 
|l.l| Fix V G Vo and let Jy be the operator on 1? given as multiplication by sgn(V) (where, 
for definiteness, we can set sgn(V)(x) = +1 whenever V(x) = 0). Then J v = Jy 1 = Jy while 



v = v\v\JvV\v\. 

Define an operator Ay by 



Ay = y/\V\T -W\V\. 



Using similar ideas to those in the proofs of Theorem |1.1| and Lemma 3.1 we get the following 



Lemma 3.2. The operator Ay is a compact self-adjoint operator on L 2 . Furthermore, we have 

7 G Xy iff — 1/7 G spec(JyAy). 

If V is single-signed we can choose Jy = ±7 (+1 if V > or —I if V < 0). Then 7 G Sy 
iff — 1/7 is in the spectrum of the compact self-adjoint operator ±Ay. This gives a justification 



of Theorem 1.5 although a more elementary argument is also possible (see Section 3.2) 



Remark 3.3. Lemma 3.2 can be viewed as a Birman-Schwinger principle for the Dirac operator 
Ty. This wide ranging principle has been used to obtain a number of results related to those 
presented here, both in the single-sign case where the associated Birman-Schwinger operator Ay 
is self-adjoint (see [Kl , IBi La] ) . and in the variable sign case where we need to consider the non- 
self-adjoint operator JyAy (see [5"a]). Approaches based on the Birman-Schwinger principle 
rely on obtaining spectral information about the operator Ay or JyAy. Potential sources of 
information include eigenvalue or singular value estimates (such as [Cw]), or pseudo-differential 



techniques leading to eigenvalue asymptotics (such as [BiSoJ). In Section 3.3 we take a different 



approach, based on Priifer techniques; this is convenient in our one-dimensional setting, and 



allows for slightly less restrictive assumptions on the potential (see Remark 3.16 in particular) 



For any V G Vq note that \V\ is a single-signed potential which is also in Vq. Using the 



association with compact operators given by Lemma 3.2 we are able to link points in Sy and 
Siyi with the eigenvalues and singular values of a single operator, and thus estimate the former 
using the latter (via Weyl's Inequality). Let < p,\ < 1x2 < ■ ■ ■ be the positive points in Siyi, 
ordered by size. Also let 71,72, • • • denote the points in Ey n {z G C : z G M + or Imz> 0} 
ordered so as to have non-decreasing modulus (and counted according to algebraic multiplicity). 
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Lemma 3.4. For any N G N we have 



N N 

Y[hn\ > Yl^n 

n=l n=l 

Proof. By symmetry we know that the points in £y are just ±j n , n = 1,2,..., while the 



points in Eiyi are just ±/i„, n = 1, 2, . . . . Lemma 3.2 now implies ±l/7 n , n = 1, 2, . . . are 



the non-zero eigenvalues of JyAy, while ±l/[i n , n = 1,2, ... are the non-zero eigenvalues of 
J|y|A|y| = Ay (note that, we can take J\ v \ = I). However (JyAy)* JyAy = AyJ v Ay = A v 
so the singular values of JyAy are just the eigenvalues of \Ay\; that is, l//i n , n = 1,2,..., 
where each eigenvalue has multiplicity 2. Given N G N we can now use Weyl's Inequality (see 
|Wey|) to compare the largest 2N eigenvalues and singular values of the compact operator 
JyAy, this gives 

'' 1 \2 4L, 1 \2 



n ft) 

n=l 



n=l 

The result follows. □ 

3.2 Symmetries 

Our unperturbed operator is 

T = -«t 2 V + ka 3 
Since V is real valued we immediately get 



k -V 
V —k 



(T + 7^)V = (T + 7^)^- (10) 

It follows that 7 G £y iff 7 G Sy. On the other hand, the commutator properties of the Pauli 
matrices (namely (TjUk = —Ck^j if j ^ &) give us 

cti(T + 7^) = -(T -7^)0-1 and ct 2 (T + jV) = (-ia 2 V - a 3 k + jV)<r 2 . 

From the first identity we get 7 G Sy iff —7 G Sy, while the second shows that £y is invariant 
if we replace k with —k in the definition of To. 

Remark 3.5. The symmetry corresponding to 03 can be used to help study even potentials 
(compare with our consideration of anti-symmetric potentials below). 

Proof of Theorem \l.b\ Suppose T y ytp = (To + ^V)tp = for some 7 G C and ?/> G H . 
Then (T ip, ip) = -^(Vif), if>) while (T i/>, tp), (Vif), ip) G M. If 7 ^ M we must therefore have 
(Vtpjtp) = 0. Since V is single-signed it follows that Vip = 0, leading to Toij> = —'yVip = 
and thus ip = (recall that To : H 1 — > L 2 is an isomorphism). □ 

Proof of Theorem \1.7\ We consider two symmetries of the operator T 7 y; define an anti-linear 
operator C and a unitary operator S on I? by 

Cy> = V> and (Sip)(x) = a 2 ij)(-x), x G R. 

These operators map i/ 1 (isometrically) onto il 1 , and satisfy C 2 = I = S 2 and C5 = —SC. 



Furthermore, (10) can be rewritten as CT 7 y = T 7 yC, while VS = -SV and SV = -VS (as 



V is anti-symmetric) which leads to 5T 7 y = — T^yS. 
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Now suppose 7 G Sy n M and choose ^ ift G i? 1 which satisfies T 7 yj/> = 0. Then 
T 7 yCV> = CT 7 yV = and T 7 yflty = -ST 7 y</> = 0. 



However is a simple eigenvalue of the operator T 7 y by Lemma pTlj so we must have Crp = aip 
and Stp = for some a, (3 G C. Then 

H 2 V = C 2 -0 = V, #V = S 2 ip = tp and a^ = = = -a/3t/>, 

so | ce | 2 = 1 = /3 2 and a/3 = —a/3 (recall that ^ ^ 0). These equations clearly have no solution, 
so we must have Y,y n E = 0. □ 

3.3 General bounds and asymptotics 

Suppose V G L\ oc . We can view our basic equation T 7 yi/> = as the 2x2 system of first 
order ordinary differential equations on E given by 

V^i = (fc- 7 T0V>2, 
V^j = (* + 7% 

The basic theory for such equations is well established (see [Ha]; here, and in subsequent 
references to [Ha], some straightforward modifications of the results are needed in order to 
cover L\ QC coefficients). In particular, if xq G E and ai,a2 G C, then there exists a unique 



absolutely continuous solution to (11) on E with ipj(xo) = aj for j = 1,2. Furthermore, for 
given x G E, this solution depends continuously on ai,a%, A and V (the latter as a function in 
L 1 (J) where / is the interval between xq and x). A consequence of the uniqueness of solutions 



is that for any non-trivial solution i/; of (11) we have ij>(x) / for all x G 



Now suppose 7 G E. It follows that all coefficients in (jllj) are real, so we may restrict our 
attention to solutions which are also real valued. Since any non-trivial solution ij) is absolutely 
continuous and satisfies ip(x) / for all x G E we can define an absolutely continuous function 

S : E -»• S 1 by 

, , = ^i(a) +1^2(3?) 

w + ^(*))Va 

(here S 1 denotes the unit circle in C). By lifting to E (the universal cover of S 1 ) we can define 
a further absolutely continuous function 9 : E — >• E so that 

e .» = s= *±* ; (12) 

this function is the Prufer argument associated to ip and is unique up to the addition of a 
constant in 2-7rZ. We note that Prufer coordinates are a standard tool for problems of this kind 
(see, for example, fSchj ). 

A straightforward calculation gives 



Using (11) it follows that 



™ ^lV^ 2 - V>2VV>1 - $f - ^2 

V0 = 7~ "Tj = 71/ + K —k ~r\ ' 
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However 



Ree 2ie = cos(2#) 



^ + V 2 

so satisfies the first order non-linear equation 

V9 = -fV + kcos(29). (13) 

Conversely, when V G L[ oc the differential equation ( [T3| ) has a unique absolutely continuous 
solution valid on E for a given value of #(xrj), x o G E (see [Ha]). If we have one solution 9 of 



(13) then 9 + nir provides a further solution for any n G Z (note that, for given by (12) the 



solution + 7r corresponds to taking — if) as a solution of (11)) 



Now suppose V G Vi and -0 is a non-trivial solution of (1 1 lh . Let 9 be given by (12) 



The fact that / K |V(a;)| dx < +oo means that, to leading order as x — > ±oo, ip behaves like a 
solution to (11) with V = 0. The corresponding asymptotic behaviour of 9 can be summarised 



as follows (see Section 4.1 for more details): 

Lemma 3.6. The quantity 9{x) has well defined limits as x — 

6>(±oo) := lim 9(x) G ^ + - 
Furthermore, upon restriction we have ^ G L 2 (M ± ) /ff 

7T 



±oo which satisfy 



(14) 



For any 7 G 
boundary conditions 



0(±oo) G +vrZ. 
we can uniquely specify two solutions # 7 + and 7) 



6» 7)± (±oo) = =F 



7T 



to (13) by imposing the 
(15) 



3.6 



shows that 7i ± correspond to solutions of (11) which are in L on M. . We get an 



Lemma 

L 2 solution on the whole of E precisely when these solutions 'match' at one (or, equivalently, 
any) point of E. Choosing as the point at which we check, the matching condition is simply 



that 7 ,+ (O) and 
setting 



.(0) must differ by a multiple of ir. Define a function 



Ay( 7 ) = (0 7 ,+ (+oo) - 7i+ (O)) + (0 7 ,_(O) - fty.-C-oo)) 

= -|-e 7 ,+(o) + e 7i _(o). 

We can thus characterise points in n E as follows: 



by 



(16) 



Proposition 3.7. We /7a\/e 7 G XV H 



7T 



/fl r Ay(7) G - + vrZ. 



We have Ay(0) = (note that 0n,± ar s constant functions) while Ay( 7 ) depends con- 
tinuously on 7 (this essentially follows from standard results for the continuous dependence on 
parameters of solutions of ordinary differential equations; see [Ha]). For large 7 we have the 



following asymptotic behaviour (the proof is given in Section 4.2) 
Proposition 3.8. IfV G Vi then 

A y (7) = 7 / V{x)dx + o(-f) 



as 7 



00. 
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Together Propositions |3.7| and |3.8| allow us to establish Theorem |1.4 



Proof of Theorem \1.4[ Let /jjCl denote the closed interval with endpoints Ay(0) = and 
Ay(i?). By the Intermediate Value Theorem Ay(7) takes each value in Ir at least once for 
some 7 G [0, J?]. From Proposition |3.7| we then get 



(17) 



#(£y n [0, R}) = #{ 7 € [0, fl] : Ay( 7 ) G (Z + 1/2)tt} 

1^1 1 



>#(/ R n(Z + l/2)vr) > 



7T 



2' 



where |Ir| is the length of Ir. On the other hand, Proposition 3.8 gives 



\I R \ = \A V (R)\=R 



V(x) dx 



o(R) 



as R 



oo. 



(18) 



□ 



Remark 3.9. For compactly supported potentials, Theorem |1.4| admits a simpler proof, similarly 
to |Schl Theorem 3]. The details of this argument have been written down by Michael Morin 
(supported by the NSERC Undegraduate Summer Research Award, 2011). 



If V € Vi is single-signed then Ay is monotonic (the proof is given in Section 3.4) 



Proposition 3.10. Suppose V G Vi is single-signed and non-trivial. Then Ay is strictly 
increasing ifV>0 and strictly decreasing ifV < 0. 



Using this monotonicity we can sharpen the argument given for Theorem 1.4 to establish 
Theorem 11.61 



Proof of Theorem \l.d If V G Vi is single-signed it follows from Proposition 3.10 that there 
is at most one 7 with Ay (7) = (n + l/2)vr for any given n G Z. Thus (17) in the proof of 



Theorem 1.4 can be modified to give 



#(Sy n [0, R]) = #{7 G [0, R] : Ay( 7 ) G (Z + 1/2)tt} 
= #(7 fl n(Z + l/2)7r) < ^ + 

7T Z 



The required upper asymptotic bound on #(XVn [0, R]) can then be obtained using (18). □ 



To give uniform estimates for Ay(7) we firstly introduce an auxiliary function. For any 
a > let [a\ denote the largest integer not exceeding a; that is, 

[a\ = max{n G Z : n < a}. 

Now define a function h : [0, +00) — > [0, +00) by 

h(a) 



IT 

a+ 2 



2 

—a 

7T . 



In particular, h is strictly increasing, h(a) = a for a < tt/2, a < h(a) < 2a for all a > 0, and 
h(a) + h(b) < h(a + b) for all a, b > (note that [a\ + [b\ < [a + b\ ). 



As a general bound on Ay (7) we have the following (the proof is given in Section 4.2) 
Proposition 3.11. IfV G Vi and 7 6 IR then |Ay( 7 )| < h(\j\ \\V\\ L i). 
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This result leads to uniform lower bounds for points in £y, and in turn helps to justify 
Theorem 1.3 We firstly deal with the case when V is single-signed (Proposition 3.12), and 
then consider arbitrary potentials (Proposition 3.13) by using Lemma 3.4 to reduce this to the 
single-sign case. 

Proposition 3.12. Suppose V G Vi is single-signed and non-trivial. Let < 71 < 72 < . . . 
denote the sequence of positive points in £y, arranged in order of increasing size. Then 



In > 



TT 



27 



n, 



n 



L 1 



1,2,.... 



Note that the constant in the bound is only half the asymptotic value (see Theorem 1.6). 

~>of. Suppose V > (the case V < can be treated similarly). By Propositions 3.7 and 3.1 
have Ay(7„) = (n — l/2)7r for n = 1, 2, Now if a < mr/2 then [2a/7rJ < n so 



h(a) < a + (n - 1)| < (n - -)tt = Ay( 7n ) < h(j n \\V\\ L i) 



by Proposition 3.11 However h is strictly increasing so we must have a < 7 n ||V||n- Taking 
a — > (rnr/2)~ now gives the result. □ 

Proposition 3.13. Suppose 7 £ Vi is non-trivial and let 71,72,... denote the points in 
Sy n {z G C : z G M + orlmz > 0} ordered so as to have non-decreasing modulus (and 
counted according to algebraic multiplicity). Then 



\ln\ > 



TT 



n, n = 1, 2, . . . . 



Proof. Let < //x < A*2 < • • • denote the positive points in Eiyi, arranged in order of increasing 
size. Noting that |V| G Vj is single-signed and non-trivial, Proposition 3.12 gives p, n > vn for 

n = l,2,. 
get 



where v = tt/(2\\V\\ l i). Using Lemma 3.4 and the ordering on the 7 n 's we now 

N N 



I7JV 



>Uhn\>U»n>V N N\>V N ^ 



n=l 



n=l 



e / 



for any iV G N. 



□ 



Proof of Theorem \1.3i . \f V = then Sy = and there is nothing to prove. Now suppose V 
is non-trivial. Let R > 0, set N = #(Sy n {z G C : \z\ < R}) and suppose N > 0. Using the 
symmetry Sy = — Sy we know that N = 2M for some M G N. With 71,72, • • . defined as in 
Proposition 3.13 it follows that \ jm\ < R- We then obtain 

2e 4e 
N = 2M < 2 — ||F|ki|7jif| < — ||7|| L ii2 

TT TT 

using this result. □ 

Remark 3.14. Various other estimates can obtained from straightforward modifications to the 
proof of Theorem L3 presented above. For example, if V G Vi is single-signed we can use 
Proposition |3.12| in place of Proposition |3. 13| to obtain the uniform upper bound 



#(Syn[o,i?]) <-\\v\\ L iR 



TT 
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for any R > 0. Alternatively, for any V € Vi we can estimate /i n in the proof of Proposition 



3.13 using Theorem 1.6 instead of Proposition 3.12 this leads to the asymptotic bound 



2c 



#(S v n{zeC: \z\ < R}) < — \\V\\ L iR + o(R) 



TT 



as R — y oo. 

As a direct application of the main result in |EITa| we can extend the general upper bound 



given by Theorem 1.3 to potentials in Vo n LP for 1 < p < oo. 
Theorem 3.15. Suppose V G Vo n LP for some 1 < p < oo. Then 

n {z E C : \z\ < R}) < C p \\V\\ P LP R P 
for any R > 0, where C p is a constant depending only on p. 



Remark 3.16. The main result in [EITaj is based on singular value estimates from [Cw]. We 
can't extend the argument to cover p = 1 as it corresponds to an excluded boundary case in 

3.4 Derivatives 

To study the derivative of Ay (7) we need to obtain more information about the 7 dependence 



of solutions to (13). Firstly consider any closed bounded interval I = [a, b] C M and potential 



V G L\ oc on /. For each 7 6 R suppose we have a solution 6> 7 of (13) where # 7 (6) depends 
differentiably on 7. Standard results for ordinary differential equations (see [Ha]) then imply 
6 1 {x) is differentiable in 7 for each x G / (in fact 6* 7 will depend analytically on 7 provided 
# 7 (6) does), so we can set 

GI7 



From (13) we immediately get 



Vw 7 = V - 2/csin(26> 7 ) u. 



7' 



or V(e G w 7 ) = e V, where G is any function satisfying VG(i) = 2k sin(20 7 (x)). Thus 
u, 7 (a) = e G(fe) ~ G(a) a; 7 (6) - / e GW " G(a V(x) 



7 w 7 (6) 



"^^(x) dx, 



where, for any interval J C J, 



sin(20 7 (x)) dx. 



(19) 
(20) 



We need to consider (19) in the limit b — > +00 when we take 9 1 = # 7 ,+ . Let w 7i+ denote 
the corresponding derivative in this case. Since # 7 + (+oo) = — 7r/4 is constant (recall (15)) we 
would expect w 7)+ (+oo) = 0. Furthermore sin(20 7i+ (x)) — > — 1 as x — > 00, so e 2m i a ^ — y 
as x — ^ 00. The precise properties that we require are given in the next result; these can be 
justified by straightforward if somewhat lengthy arguments. 
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Proposition 3.17. The solution # 7i+ to (13) depends on 7 differentiably while the derivative 
satisfies 

d 



w 7 ,+ (a) := — 7 ,+ (a) 



for a// a G 



There is a corresponding result for 9 7i -. Proposition 3.10 is now an easy corollary of these 
results. 



Proof of Proposition 3.10 From (16) we get 

d 



d7 



A y ( 7 ) = -w 7i+ (0)+w 7i _(0). 



Now suppose V > (the case V < can be handled similarly). By Proposition 3.17 we have 



-w 7 ,+ (0) 



i°^V(x) dx. 



Since e 2k ^ J > for any (bounded) interval JC1 the right hand side is non-negative and equal 
to only if V = on M + (as an L 1 function). A similar argument shows that w 7 _(0) is also 
non-negative and equal to only if V = on IR - . The result follows. □ 



3.5 No gaps 

For a potential V G L\ oc and interval I C E let varr(V) denote the total variation of V on I. 
We also say that V has no gaps on I if 

|i"nF -1 (o)| =0. 

To work with no-gap potentials we need estimates for the integrals of cos(6 l 7 (x)) and 
sin(# 7 (x)); to complement (see (20)) set 



<]>,/ = J cos(20 7 (x)) dx 



for any interval JCR. For large 7 equation (13) suggests 9 1 {x) should be changing rapidly 



wherever V(x) ^ 0; it follows that cos(0 7 (x)) and sin(# 7 (x)) should be rapidly oscillating, 
leading to cancellation in the integrals defining <I>j and \P j. This idea lies at the heart of the 



following result (the proof is given in Section 4.3) 



Proposition 3.18. Let I C M. be a closed and bounded interval, and suppose a potential V 
satisfies vavi(V) < +00 and has no gaps on I. Also suppose # 7 satisfies (13) on I. For any 



sub-interval J C I we have <£j, = o(l) as 7 — > 00, uniformly in J and possible choices of 
(the initial condition for) the solution # 7 . 

These estimates for $j and lead directly to the following asymptotic information about 
the change in the value of # 7 across I: 



Proposition 3.19. Let I, V and 6> 7 be as in Proposition 3.18 Write I = [a,b]. Then 



(b) - 7 (a) = 7 / V(x) dx + o(l) 



(21) 
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as 7 — > oo. Furthermore, if d6* 7 (&)/ d7 exists and is bounded in 7 then 



(0 7 (6) - 7 (a)) = / V{x) dx + o(l) 



d7 



as 7 — )• 00; /n particular, d# 7 (a)/ (I7 a/so ex/'sts and /'s bounded in 7. 



(22) 



Proof. Integrating (13) gives 



V8 1 (x)dx = >y V(x)dx + 



so (21) follows directly from Proposition 3.18 On the other hand, (19) leads to 



_d_ 

d.7 



(0 7 (6)-0 7 (a)) 



However Proposition 



3.18 



V{x) dx = u)-y(b) — w 7 (a) — / V(x)d:c 

J a 

= (l-e> 7 (6)- Al-< 

J a 

gives us = o(l) and hence 1 — e 2fc * J 



2fc* 



for all sub-intervals J C I; (22) follows. 



o(l) as 7 — ^ 00, uniformly 

□ 



Suppose a potential V 6 has support contained in the bounded interval I = [a, 6]. 
Clearly constant functions taking values in 7r/4 + 7rZ/2 solve (13) outside /. The condition 
(15) (together with the uniqueness and continuity of # 7 ,±) then gives us 6 y -(x) = 7r/4 for 
x < a and 6> 7i+ (x) = — 7r/4 for x > b (this can also be seen from the form of the corresponding 
solutions to (11)). When defining Ay in this case it is convenient to choose the left endpoint 
of I as the point at which to check whether # 7 _ and 6* 7j+ can be 'matched'; alternatively, we 
can keep our current definition of Av if we simply translate our problem so that a = (see 
Remark 1.8). Making such a choice we get 

Ay (7) 



7T 

2 



,+(a) + 7 _(a) = 9 7:+ (b) - 7) +(a), 



where 7j+ (6) = — 7r/4 is constant (as a function of 7). We will work with Ay in this form for 
the remainder of the present section. 

For compactly supported potentials without gaps we can now rephrase the conclusions of 
Proposition |3.19 to get the following improved and extended version of Proposition 3.8 



Corollary 3.20. Suppose V G BVq has no gaps. Then 

Ay(7) = 7 / V(x) dx + o(l) and Ay( 7 ) = I V{x) dx + o(l) 

JR d 7 JR 



as 7 



00. 



When f R V(x) dx 7^ it follows that Ay is monotonic for sufficiently large 7. In this case 
Theorem |1 .9| can be proved using an argument very similar to that used for Theorem 1.6 The 
case J R V(x)dx = can be treated with a separate observation. 



3.20 



enables us to choose S > so that 



Proof of Theorem \l.G{ If f R V(x)dx = Corollary 
I Ay (7) I < 7r/2 whenever \ j\ > S. Then Sy n M. C [-S, S] by Proposition 3.7 However Sy is 
a discrete subset of C (Theorem |1.1[ ) so Sy n [—S,S] contains at most finitely many points. 
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Now suppose f R V(x)dx / 0. By Corollary 3.20 there exists S > such that Ay(7) is 
strictly monotonic for all 7 > S. Suppose R > S and let Is,r denote the closed interval with 

O(l). 



endpoints Ay(£) and Ay(i?). Arguing as for Theorem 1.6 we then get 

\Is,r\ 



However Corollary 3.20 



n [S, R}) = #(l S)R n(z + i/2)7r) 

also gives 



7T 



\I s ,r\ = \A V (R) - A V (S)\ = \A V (R)\ + O(l) = i? 



+ 0(1) 



as i? — > 00 (note that, 5" is fixed). The fact that £y n [0,5] contains at most finitely many 
points (see above) completes the argument. □ 



4 Technicalities 



4.1 Asymptotic behaviour of ode solutions 

We return to viewing our basic equation T 7 yi/> = as the system of ordinary differential 
equations (11). When V = we have the exponential solutions 



1 

±1 



(23) 



(see Lemma 2.1). When V G Vo (or, more generally, V belongs to the locally L 1 version of Vo) 
we can find solutions of (11) with similar asymptotic properties asx-> ±00 (see [Ha], Chapter 



X). In particular, there are non-trivial solutions ip± and ip± which satisfy 



lim x 

x— ¥+00 



1 log|y>j|(a;)| = ±k and 



lim x 

x— y— 00 



_1 log|(£>_j_(x) I = ±k. 



Since k / tpX( x ) anc ' f-( x ) nave different asymptotic behaviour as x — > +00 (the growth 
of (f±(x) as x — > +00 is roughly like e ±kx ); it follows that these solutions must be linearly 
independent. A similar discussion applies to <p+ and (pZ- 

Lemma 4.1. IfV G Vq then either ^ spec(Ty) or is a simple eigenvalue ofTy. 



Proof. We have G spec(Ty) iff is an isolated eigenvalue of Ty (see Remark 1.2). Now 



suppose ip is an eigenfunction corresponding to 0. Then if> satisfies ( 11 ) (with 7 = 1). Since ip+ 
and <p! are linearly independent solutions of this equation we must have ip = a+ip^+ct-ipZ for 
some constants a±. Restricting to the interval M + we have <£>t ^ L 2 (IR + ) while ipZ G L 2 (M + ). 
It follows that a+ = 0, and so ip is a multiple of ip_. This must also be true for any other 
eigenfunction corresponding to 0, so any two such eigenfunctions are linearly dependent. □ 

Remark 4.2. Using a similar argument we can also get = /3+ipZ for some constant 
showing that ipZ and <p~Z are linearly dependent. In fact this is an alternative characterisation 
of when is an eigenvalue of Ty. 

When V G Vi solutions to (11) have well defined leading order asymptotics as x — > ±00; 



these asymptotics are solutions to the same equation with V = 0, so must be linear combinations 



of the exponential functions given in (23) (see [Ha], Chapter X). In particular, we can choose 
(p± and (p± so that 

' 1 



lim 



>Tkx 



-Aikx 



±1 







(24a) 
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while 



lim 



^nkx 



±1 



0. 



(24b) 



Furthermore ipZ and ipZ are uniquely determined by these asymptotic conditions. (The solutions 
tp\_ and ipZ are however only determined up to the addition of a multiple of <pZ and (fZ 
respectively.) 

Proof of Lemma LT31 We consider the case x — > +oo; x — > — oo can be handled similarly. Since 
and cpZ are linearly independent we can write = a + ip+ + a-ipZ for some constants a± 
(which can't both be 0). Now by ([24a]) 



,„.., (a +¥3 +(x) + a_^(x))i Jl ifa+^0 

lim = lim = \ 

x^+oo rj) 2 (x) %-*+oa (a+(pT(x) + a-<pZ{x))2 1—1 if a+ = 0. 



tpi(x) 



However tan6> = ipi/ipz so ([Mj) (for s -> +oo) follows. On the other hand, ip G L 2 (M + ) iff 
a + = 0, leading to the second part of the result. □ 



4.2 Estimates and asymptotics for A t 



Propositions |3.11| and |3.8| both follow almost directly from the following estimate: 

—it/ '4. For any x G E we 



Lemma 4.3. /.et V G Vi anc/ suppose 9 solves (13) w/t/7 #(+oo 

(x) + ^|<^| 7 |^ O °|y ( t)|di 

4 similar result holds on (— oo, x] i/i//?e/i 0(— oo) = 7r/4. 

The basic idea is that if Q(x) increases on an interval [a, b] crossing a range of values where 
kcos(29(x)) is non-positive then we must have 



jV(x) dx > 9(b) -9(a). 



We can then add the contributions from each such crossing. 

Proof. Let n G Z and suppose 6>(x) G [—(n + 3/4)7r, — (n + l/4)7r] for all x in some interval 
[a, b] C M. Then fc cos(20(a?)) < so V<9(x) < 7^(1), and hence 



9(b) -9(a) < / jV(t)dt 



(25) 



This estimate continues to hold with the obvious interpretation when b = +00. 

Now let x G R and suppose 9(x) < —it/A. Choose n G No so that 6(x) G (— (n + 
5/4)7r, — (n + l/4)7r] . Using the continuity of 9, and the assumption that 6>(+oo) = — 7r/4, we 
can now choose a sequence of points 

x < a n <b n < a n _i < 6 n _i < • • • < a < 6 < +00 

such that 

(i) For j = 0, . . . ,n and x G [a 3 -,6 3 -] we have 9(x) G [-(j + 3/4) /7r, -(j + l/4)7r] . 
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(ii) For j = 0,...,n we have 9(a,j) = -(j + 3/4)ir and 0(bj) = -(j + l/4)7r, with the 
exception that 9{a n ) = 8{x) if 9{x) G (-(n + 3/4)7r, -(n + 1/4)tt] . 



Applying (25) we then get 

11 n rbj roc 

E( e ( 6 i)- e ( a i))<E/ 7^)dt</ i-yn*)idt. 



(26) 



j=0 



j=0 " a J 



If 0(x) G (—(n + 5/4)7r, — (n + 3/4)tt1 then each term in the sum on the left of (26) is ir/2 so 



7T 



(n + 1)- < / | T T^(t)|d*. 



It follows that 



(n + 1)- < - 
v '2 ~ 2 



l7^(t)|dt 



On the other hand 9{x) > — (n + 5/4)tt so 



7T 



— — <(«+ l)7r < h ( / |7V(t)| dt 



Alternatively suppose 8(x) G (— (n + 3/4)7r, — (n + 1/4) 7r] so 9{a n ) = 9{x). Then the term in 
the sum on the left of (26) is tt/2 for j = 0, . . . , n — 1, so 



n 



n 1 

2=E(%)-%))< / 17^(4)1 dt, 



leading to 



Furthermore 



3=0 



7T 7T 

n— < — 
2 ~ 2 



hv(t)\dt 



J>(6,) - 9{ aj )) = - (n + -) tt - 9{a n ) + n| = -0(x) 



3=0 



7T 7T 

n—. 

4 2 



Using (26) again gives 



7T 



0(x) " T < / hV(t)\ dt+n-<h( \ 7 V(t)\ dt 



A similar argument can be used to deal with the case 9{x) > — 7r/4 (we need to consider 
intervals where 9{x) decreases across the range [(n— l/4)7r, (n + l/4)7r]). □ 



Proof of Proposition \3.11 By Lemma 4.3 we get 



7T 



fty,+(0) + T <h h\ \V(t)\dt . 



The equivalent estimate on (— oo,0] gives 



V(o)-||<MH / |V(/)|d/ 



The result now follows from (16) and the fact that h(a) + h(b) < h(a + b) for any a, b > 0. □ 
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Proof of Proposition 3J3 Recalling (16) we can write 

A y ( 7 ) = A+( 7 ) + A y ( 7 ). 



where Ay (7) := — 7r/4 =f ^ 7 ,± (0). Now suppose K > 0. Since # 7i+ satisfies (13) we get 

r-K t-K 

6 J>+ (K) - 6> 7l+ (0) = 7 / V(x)dx + k cos(2(9 7i+ (x)) dx 



so 



On the other hand Lemma 4.3 



7 , + (iT) -0 7 ,+ (O) -7 / F(x)dx 
gives us 



< | 7 | / |V(x)|dx + jfeif. 



7T 



7)+ (if) + j <h[\-y\ I \V(x)\ dx ) < 2| 7 | I \V(x)\ dx. 



K 



K 



Therefore 



A+(7)-7/ V(x)dx 



e i!+ {K) -J + 7j+ (tf) - 7 ,+ (O) - 7^" F(x) dx 

/■oo 

<3| 7 | / |y(x)|dx + fcET. 



Now choose if = K 7 > for each 7 so that K 7 — >• 00 and 



as 



00 (we can 



take K 1 ~ 17^ with < p, < 1, for example). Then J^ |V(x)| dx = o(l) (since V G L 1 ) so 

/>oo 

3| 7 | / |F(x)|dx + fcK 7 = 0(7) 
— > 00. A similar estimate can be obtained for Ay( 7 ). □ 



as 



4.3 Estimates for $j and \Pj 
Intervals 

It is easiest to deal with the potential V in pieces where it is single-signed and bounded away 
from 0; the next result is the key to identifying these pieces and sets up some of our notation. 

Lemma 4.4. Suppose V G BVq has no gaps and set I = supp(V). For each e > there exists 
a finite collection of disjoint closed intervals J e n C I, n G M e , such that 

(i) For each n G N e V has constant sign on I Ei7l and \V(x)\ > e for all x G I e>n . 

(ii) Setting E £ = I \ U„ e AT £ T e,n we have \E £ \ -> as e -> 0. 
(Hi) #A4 < ve^ 1 , where v := vav^(V) is the total variation ofV. 

We would ideally like E £ to be V~ 1 {{—e,e)), so the union of the I £l n's would be 

{x : V(x) < -e}U{x : V(x) > s}. (27) 
There are however several technical issues associated with this choice; 
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• We are not assuming that V is continuous so (27) may not be closed (or even particularly 
well behaved). 

• Even if V is continuous the number of intervals in V~ l {(— e,e)) may be uncontrollable 
(even infinite); we need a useful bound on this quantity. 

The technicalities in the following argument arise from the need to deal with these issues. 
Proof of Lemma\4jl Set 

Xf = [x — e, x + e) and X~ = (x — e,x + e); 

xeV^d-oo^e)) xGV' 1 ((-2e,+oo)) 

(so Xf is the 'e-neighbourhood' of {x : V(x) < 2e}, while X~ is the 'e-neighbourhood' of 
{x : V(x) > -2e}). 

Since X+ is an open subset of R it consists of a countable union of disjoint open intervals; 
let F £ + C be the union of those intervals which intersect y _1 ((— oo, e)). If J = (a,/3) is 
a maximal interval in Y e + it follows that we can find x G J with F(x) < e. Furthermore J 
must also be a maximal interval in so, if J is bounded, its endpoints satisfy a,/3 ^ X+. 
Thus V(a), VC9) > 2e and so var 7 (V) > \V{x) - V(a)\ + \V(fi) - V{x)\ > 2e. If J is a 
half-infinite interval a similar argument shows that varj(V^) > e. 

Now let M+ and iV+ denote the number of semi-infinite and bounded maximal intervals in 
Y+ (we're presently allowing N+ = +00; if Y+ = R set M+ = 2, JV+ = -1). Then 

v = var M (y) > var-+ (V) > eM+ + 2eiV+ => N+ < 1 (z/e -1 - M+); 

ie 2 

in particular, must be finite. However V has compact support, so Y e + must be unbounded 
both above and below. Since Y+ consists of a finite collection of intervals, it must therefore 
contain semi-infinite intervals at either end; that is, M+ = 2. The set Zf := R \ Y e + is then 
the union of iV+ + 1 closed bounded intervals; write Zf = UneA/ 4 ' ^ £ > n wnere ' s some 
indexing set with #A/" e + = + 1. It is straightforward to check that we have V(x) > e for 
any x G Zf. 

We can similarly define Y~ , Z~ , I £ , n for n G Af~ and AT". Set M £ = 7V+ UJV". Property 
(i) is immediate, while (iii) holds since 

#A4 = #A4+ + #N~ = N+ + N~ + 2 < 2^e- 1 - 2) + 2 = i/e" 1 . 

Now UneA/i ^ £ > n = U ^e" =: If < £ i < £ 2 a straightforward check gives us 

X% C X± => F± C F± => Z± D Z±, 

so Z £l 5 Z £2 and hence £ , £l = I \ Z £l C / \ Z £2 = E^. Since J5 e C J for all e > and 
|/| < +oo, property (ii) will now follow if we can show |P| e>0 .E £ | = (see [Ru], for example). 

For each 5 > let 

^ = {^n £ >o^ = ±V(x)>5}. 

Suppose xi,. . . ,xn are distinct points in and let < e < mmi<ij<jy\xi — xj\. From the 
definition of Sl^" we can find Xj j£ with V(xj )£ ) < 2e and \xj j£ — Xj\ < e for j = 1, . . . , N. It 
follows that 

JV 

v = var R (V) > Yy( x i,e) ~ V ( x i)\ ^ N ( S ~ 2£ )- 
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Taking e — > + gives N5 < v. Hence Qf is finite (with #£lf < v5 x ). Clearly we also have 
n+ DO+ if < <5i <5 2 . 

Similar properties hold for Qj. It follows that the set 



O:=nWnx-)\y- 1 (0)c(J(o+u^) = \J(n+ n un- /n ) 

e>0 8>0 nSN 

is countable (it is contained in a countable union of finite sets). However 

E e = i\z £ = i\(z+uz-)QY+n y- c x+ n x - 

so 

n ^ c / n p| (x+ nl £ -)c(/n t^ro)) u n. 

e>0 e>0 

The no gap condition on V and the countability of $7 now imply |P| £>0 .E £ | =0. □ 
Integrals 



To justify Proposition |3. 18| we start by considering the cancellations over each period of cos or 
sin (Lemma 4.5) and deal with any incomplete periods (Lemma 4.6). In both cases we work on 
an interval J = [arj,ai] where V9{x) > 0, so 9 is invertible. Making the substitution u = 29 (x) 
we get 

dn = 2V9(x) dx = 2[jV(x) + kcos(29(x))} dx, 

so 

r a l f29(a 1 ) f( \ 

f(29(x))dx= / r — r . ,\ — dn (28) 

7 2e(ao) 2{ 1 V{9-\u/2)) + k C o S {u)\ K ' 



for any / (we will take either / = sin or / = cos). 



Lemma 4.5. Suppose 9 satisfies (13) on an interval J = [acb a i] where V(x) > e with 7 e > k. 
If9(a 1 ) = 9(a ) + n then 

k\J\ vaij(V) , IT . varj(y) 

< , 1 1 + , \\ and \Vj\< 



2(76 — k) e(7£ — k) e(-)E — k) ' 

Note that we get a better estimate for the sin integral; the extra term in the estimate for 
the cos integral is needed to cope with the fact that this integral is non-zero even when V is 
constant. 

Proof. Since 9{a\) = 9(ao) + it we can define x '■ [0, 27r) — > [29(ao), 20{a\)) to be the unique 
bijection with x( u ) ~ u S 2-7rZ for all u G [0, 27r); \ ' s piecewise affine with at most one jump 
(x will have no jumps iff 9(o,q) G 27rZ). Using ([28]) we can then write 



ai f 27T sinful 



fyj = / sin(26>(x)) dx = I r t^a-i ? ? \ /o(\ 1 ? ~ dn = ^+ + 

J a 



where 

r sin ( n ) 



o 2[ 7 F(e- 1 (x(«)/2)) + fccos(u)] 

r sin(n) 

o 2[ 1 V(^(u)) + kcos(u)] 



with £+(n) = 9- 1 (x(u)/2), and 

sin(v 

lo 2[ 1 V(9- 1 (x(u + tt)/2)) +fccos(u + 7r)] " 7 2[ 7 F(£_(«)) - fecos(u)] 



= r sm( ^ } : ^ du = - r sin(u) . ,„ dn 
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with = 9 1 (x(w + 7r)/2); £ + and are bounded piecewise continuous functions [0, ir] — > 

M with ranges contained in J = [ao, a\\ (in fact the range of £± is just the set of x G J where 

±sin(20(x)) > 0). 

We now seek bounds on ty± using (constant) bounds on V. Set m = val{V{x) : x G J} 
and M = sup{V(x) :i£j}so6<m<M< +oo and m < V(x) < M for all x G J. Thus 

< 7m ± Acos(u) < 7V(£±(it)) ± Acos(u) < 7M ± A;cos(u), u G [0, 7r], 

so 

Z' 71 ' sin(u) , /"^ sin(it) 

7o 2[7M± fccos(n)] ~ Jo 2[ym ± k cos(u)} 

The integrals appearing in these bounds can be calculated explicitly, leading to 

1 , /7M + AA 1 , f^m + k\ 

— log — < < — log — . 

2k & \jM-kJ- ± ~2k *\Tm-k) 



Hence 



|*j| = |*+-(-*_)| 
1 



< 
~ 2k 



, /Sm + AA , /7M + A 
log - log 1 



\7m — A; / \7ikf — fc 



1 / 7(M - m) 

log 1 + 2A- 



2A; °V (jm - k)(-fM + k) 

< 7(M-m) 
~ (7m — k){^M + A;) 

(note that, log(l + i) < i for any i > 0). Since (7m — k)(^fM + k) > (je — k)je > and 
M — m < vaij(V) the required estimate for now follows. 

In a similar manner we can write 



$j = / cos(20(x)) dx = + $_ 

J an 



where 

*. = ± r 72 ^ d « 

± i_ 7r/2 2[7y(r /± (n))±Acos(n)] 

for some bounded piecewise continuous functions rj± : [— 7r/2,7r/2] — > R whose ranges are 
contained in J = [ao,ai] (the range of rj± is just the set of x G J where ±cos(26(x)) > 0). 
For u G [— 7r/2,7r/2] we have < cos(u) and 

< 7m — Acos(u) < 7V(r/±(u)) ± fccos(u) < 7M + A;cos(u), 

so 

^Z 2 cos(u) , x r' 2 cos(n) 



y_ 7r/2 2 7 M+A ; cos« - ± -y_^ 



/ 2 2[7M + Acos(n)] 7-^/2 2[7m — k cos(tt)] 

Now, for any < c < 1, (7m - kc){^M + fee) > (7m - AX7M + k) > (recall that M > m 
while 7m > 75 > k)\ hence 

Q< 1 1 2kc + ~f(M-m) < 2kc + j(M-m) 



— kc 7M + Acc (7m — Ac) (7M + Ac) (7m — k){^M + A;) 



28 



Putting these estimates together gives 



< 



< 



tt/2 



cos(u) 



_ 7r/2 2 [7m - fccos(w)] 



du 



tt/2 



cos(u) 



du 



k 



(jm-k){jM + k) J_„ /2 



tt/2 



_ w/2 2[ 7 M + kcos(u)j 

7(M - m) 



+ 



cos (u) du + 
7(M - m) 



2{ 1 m-k){ 1 M + k) J_ n/2 



tt/2 



cos (it) du 



2(7771 — k)(yM + fc) (7m - k)(yM + k) ' 



(30) 



The second term can be estimated as for (29). On the other hand, (13) implies |V0| < 7M + k 
on J, so 



vr = 9(ax) - 0(ao) < (7M + fc)(oi - a ) = (7M + k)\J\. 



The required estimate for the first term in (30) follows. 



□ 



Lemma 4.6. Suppose 9 satisfies (13) on an interval J = [acb a i] where V(x) > e with ye > k. 

lf\0(ax) -0(ao)\ <vr then 



ye — k 



Proof. For any u G [20(an), 20(ai)] we have 

7y(0 _1 (n/2)) + fc cos(tj) > ye - fc > 



(note that 9 1 {u/2) G J). For either / = sin or / = cos (28) now leads to the estimate 



f(29(x))dx 



< 



20(oi) 



du < 



29(a ) 2 (7^ - fc ) 2 (7£ - k) J 2 e(a ) 



2e(a )+2w 



\f(u)\du 



ye — k ' 



where the middle step follows since 9(a%) < 0(an) + 7r. 



□ 



The previous Lemmas can be combined in a straightforward way to deal with more general 
intervals: 



Lemma 4.7. Suppose 9 satisfies (13) on an interval J = [a,b] where V(x) > e with ye > k. 
Then 

A similar result holds in the case that V{x) < —e for all x G J. 



Proof. Since is absolutely continuous (as is any solution of ( 13 ) when V G £i oc ) and W(x) > 
76 — > we can choose n > and points a = ao < a\ < • • • < a n < b such that 
9(a,j) = 0(oj_i) +7r for j = 1, . . . ,n and 0(6) < 6(a n )+ir. Set Jj = [cij_i, oy] for j = 1, . . . , n 
and J = [a n , b] so Ji, . . . , J n , J have disjoint interiors while J = 3\ U • • • U J n U J. Now 

i*,i<£i*, j i + i* jH <£^ + 2 



^ e(7e - k) ye 



k 



using Lemmas 4.5 and 4.6 The ^ j estimate now follows from the fact that varj 1 (y) + • • • + 
Y3xj n {V) < varj(F). A similar argument leads to the $j estimate. □ 
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We can now combine the previous result with Lemma 4A to deal with arbitrary sub-intervals 

of I = supp(F): 

Proposition 4.8. Suppose V G BVq has no gaps. Let e > and consider the notation of 
Lemma 4.4 If 9 satisfies (13) on I and 76 > k then for any sub-interval J C I we have 

, Qve~ 1 + k\J\ 1 1 1 3^e _1 

Proofs Set J £ , n = I £) „ n J for n G M £ , and W £ = {n G W £ : J £ , n ^ 0}. In particular 

#K < #A4 < ^ (3i) 



by Lemma [4~4|jii), while 



F £ := J \ |J J £ ,„C/\ J I £>n = £ £ 

nSA/"' nSA/i 



so 



1^1 < \E E \. 



(32) 



Now 



sin(20(a;)) dx. 



Using ( [32]) the modulus of the final term is clearly bounded by \E e \. On the other hand, by 
Lemma [474] ^ i ) V has constant sign on J £>n C I en and satisfies |V(a:)| > e for all x G J £>n . 
Lemma |4.7| thus gives 

2 + e~ 1 var J (V) 

However the J e . n for n G A/" £ are disjoint subintervals of J C 7, so 

^ var Jen (F) < varj(F) < var K (y) = 1/. 



Together with (31) we then get 



(2 + e- 1 var Jen (V)) <2#M' e + v£~ 1 < Zve' 1 . 

neAfl 



The estimate for follows. Since 



a similar argument leads to the estimate for 



<\J\ 



□ 



Proof of Proposition \3.18\ Choose e = e 7 > for each 7 > so that e 7 — >• and e 7 7 — > 00 
as 7 — > 00 (we can take e 7 ~ 7 _A1 with < fi < 1/2, for example). Letting 7 — » 00 it follows 
that ve- 1 ^ - k)- 1 -> and |£L I (f rom Lemma |4.4[ii)); both limits are independent 

at a 

rate that can be bounded uniformly in J. The result now follows directly from the estimates in 
Proposition pT8| □ 



of J. On the other hand, |J| < |/| (since J C I) so we also have |J|(7£ 7 — k) 



-1 
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